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The spectral density of a frequency modulated carrier is evaluated for the 
case when the modulating baseband wave is a "quantized" random facsimile 
signal. By this nonsynchranous form of modulation we mean holding one 
of the allowed set of transmitted frequencies for a finite, but randomly dis- 
tributed, lime before switching to another frequency while maintaining phase 
continuity. Emphasis is given to the Poisson case of exponentially dis- 
tributed intervals between transitions, and some typical curves for discrete 
level and continuous level situations are included. 

I. INTRODUCTION 

In facsimile data transmission systems, printed or pictorial informa- 
tion is converted into electrical signals by optical means. At any instant 
the signal corresponds to a definite grey level of the facsimile copy. The 
resulting electrical wave is an analog signal and can be transmitted as 
such. In some applications only black-and-white images need be trans- 
mitted and therefore the electrical signal may be quantized into only 
two levels. If more detail is desired, multilevel quantization can be 
applied. 

It is possible to model such a quantized signal by considering a random 
sequence of points on the time (0 axis. At each point a transition may 
occur in the signal. The value of the signal between transitions is a 
constant, taking on one of N different values. For the black-and-white 
case, there are only two permissible values, either +1 or —1. The quan- 
tized facsimile signal is statistically characterized by specifying the 
distribution of the points on the t axis and the distribution of the ampli- 
tudes between transitions. 

In this paper we concern ourselves with the spectral density of a car- 
rier wave whose frequency is modulated by quantized facsimile signals. 
The spectral density is a useful item in the statistical description of such 
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a signal in that it furnishes an estimate of bandwidth requirements. It 
often is also used to evaluate mutual interference between channels. 

So far as is known, the amplitude modulation case is the only one 
hitherto covered in literature. However, we were prompted to examine 
the FM case as FM is currently used in facsimile data sets. From a 
practical point of view the most interesting case is that in which the 
phase is continuous at the transitions, as may be obtained from keying a 
single oscillator. This case differs from previous results 1 in that the 
transitions occur at random times. 

The present paper gives a complete solution for the spectrum for an 
arbitrary distribution of the interval between transitions as well as 
arbitrary distribution of amplitudes. We treat an important special case 
of Poisson transitions, for which we present our results graphically in 
terms of the important parameters of the process. 

An interesting feature is the rapidity with which the spectral density 
falls off with frequency measured from midband as compared with the 
AM case. The extent to which spectral peaking occurs at the average 
signaling frequency for some range of parameters is another curious 
feature. As would be expected from the asynchronous nature of the 
modulation there can be no steady sine-wave components in the process 
and therefore there are no discrete components in the spectrum. 

II. ANALYSIS 

The baseband facsimile signal is constructed in the following form. 
Pick a finite set { t j at random and arrange the points such that 

= to < k < h < ••• t N = T. (1) 

Define a set of functions 



/l, U £ t g tn 
where 



9K(t- <«) = j , elsewhere (2) 



A„ — t n +l — tn • 

In terms of (1) and (2), construct the baseband signal x(t) as the fol- 
lowing time series 

X(t) =£ Cfl&.C* - tn) (3) 

n=0 

where a = (ao, fli , ■ ■ • , a.v-i) is an additional arbitrary set of iden- 
tically distributed random variables. 
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The instantaneous phase \p(t) is represented as follows 
*(0 = w c t + u d [ .v(t')dl' + <p = yPx(t) + <p, 0£t£T (4) 

Jo 

where ip is uniformly distributed on [0,2r], giving the value of the phase 
at i = 0. The instantaneous frequency is d\p(t)/dt = w a + av.r (/), 
where co f is the angular carrier frequency and ov is the minimum angular 
frequency deviation. 

The FM wave whose spectrum we wish to examine has the following 
representation 

S(t) = A cos [*(*)] 

(5) 
= (il/2)exp{/<K0l + (4/2) exp {-#(01, 

where A is a real amplitude. The spectral density of S(t), G(u), defines 
the average power in a unit bandwidth. Formally, it may be obtained as 

G(u) = lim (2/T) ( | S(<a,T) | 2 >, a > 0, 
where »S'(co,7') is the Fourier transform of 5(0 given by 

S(u,T) = [ S(t) exp (-iu>t) dt (7) 

• 

and the symbol ( • ) denotes the ensemble average over all the random 
variables in S(u,T). 

One may write S(o),T) as 

S(o,T) = U/2)c ,v lF 1 (c,r) + (A/2)e-*Wi(u,T) (8) 

where FPj and IF 2 arc the Fourier transforms of exp [#1(0] and exp 
[— i\f/\(t)] respectively. 

The ensemble average of | S(u,T) |" over <p is readily carried out to 
obtain 

< | S( W ,T) | 8 >, = (A 2 /4) | W 1 (u>,T) | 2 + (A 2 /4) \ W,(co,T) | 2 . (9) 

We proceed to evaluate ( | \\\(u,T) |"'): 

Wi(a>,T) - / expzj^i(/) - ut]dt 

«-*-. ,.,„ < 10 > 

= L / exp t{^i(0 - ut\dt. 

h=n J l l: 

We observe from (4) that in the interval [t k , fc +1 ], 
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n=k-l 



$l(t) = U>d 2 a « A » + <*df*k(t — fc) + U c t. (11) 

Inserting this expression into (10), and using 

tk = 2 A »' 

n = 

we find that 

iFiwn = 5 ^ [ ex p {* £ x » a «} - ex p {* jic x » a » }] < i2 ) 

where 

X„ = LOtldn — CO + tOc ■ 

Multiplying (12) by its complex conjugate we obtain 

| W l ( u ,T) f = 2Re [*£* 1 ~ "P tXA 

L fc-o A «= 

+ E i (expjt Z X„A„i + exp(tEx„A„l (13) 

A,8=0 AA-Ag \ [ ;i=8+l J { «=« J 



A,8 = (l 



i 2 X„A„> - exp U 2 X„A„>j , 



— exp < i 

where Re ( • ) denotes the real part. 

At this point we must specify in more statistical detail the sets of 
random vectors A = (Ai , A 2 , • • • , Aw) and 7. = (Xi , X2 , • • • , Xw). In 
our original representation (3), the A„'s are the intervals between transi- 
tions. We adopt the reasonable assumption that these intervals are in- 
dependent. On the other hand the random variables X„ which are re- 
lated to the amplitude a n by (12) are not independent if we consider 
only observable transitions. Clearly for observable transitions one in- 
sists that a n 9± a n -i , and thus adjacent amplitudes are dependent. To 
remedy this awkwardness in the analysis, we construct an alternate 
random process with independent a„'s by admitting virtual transitions. 
That is, we do not require that the signal change at every t„ given in the 
sequence (1). We will show later how such a process, entirely equivalent 
to the original, may be constructed. For the present, we shall merely 
assume that the X„'s are independent. 

For fixed N the average of (13) with respect to A and 7. becomes 
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(| Wi(u t T) | 2 >aa = 2Re 



v( 'l\ N Asrt 



k,8=N—\ 

+ z 

A", 8=0 
k>a 



exp i 2 X»A„ 



\ 



\ 
exp i zL XnA« 



\ 



+ 



14) 



\ 



X.X* 



/ exp i 2 X„A n \ 

\ — / 

\ X.X* / 

/ • £ 

/ exp i 2-, A„A„ 



\ X,X* / j j 

The respective averages in (14) may be expressed in terms of the char- 
acteristic function of A. A typical calculation yields 

/ * \ 

/ exp i Yli XAn \ / 1 TT ... A V 

_ / 1 ' n=A - 



\ 



A S X*- 



3l,A 



(^ II expiX.A,) \ 

\A.sAt \, i= 8 +i / A/ }. 



(15) 



= <^V ff c A (x„)\ 

\A^Aa- n=«+l / 



_/l\ /C A (X)\ 



X/x \ X 



[<C A (X)) X ]*— ' 



where C* (X) = (exp (/XA))a is the characteristic function of the random 
variable A. 

Using the same procedure as above on every term in (14) we obtain 



(| Wiiuf) | 2 )aa = 2Re 



K^# 



i - Ca(x: 



k>a 



(16) 
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where 

P- (C A (X)>x. 

Since \ p\ < 1, the series in (16) can be summed. If we designate the 
average of N by N, divide (16) by T and take the limit as T — > °o such 
that 

v = lim (N/T) 



we obtain 
1 

T— » J 



lim ^ <| Wi(a,T) | 2 > AA,iV 



<=,Ke[L-^\-«^^xy r ^], 

where we made use of the identity 

N-l N-l 



(17) 



n-l 



k,n=0 
k>, 



We can repeat the identical operations on Wi(a,T) in (9) that we 
have just concluded on Wi(u,T) and obtain an identical expression 
except that co — w c in Wi will be replaced by w + u c . 

Combining (17) with (9) and (6), we write down the positive image 
spectrum as our general result, namely 



G+(u>) --= Ap 
Re 



/f l - Ca(X)\ 1 // 1-C,(X) \\ 2 ] 

\L x 2 /x i - (Ca(x))x \\ x /J y 



is) 



This result is general and applies when the choice of amplitudes is 
made independently at every t n point. As remarked earlier, in a real 
facsimile process the choice of amplitudes is constrained. If a transition 
is to occur at every t n point, the adjacent amplitudes must be correlated. 
We now show that the real process can indeed be represented in terms of 
independent amplitudes by the expediency of introducing virtual transi- 
tions. We write down the following equality 

E ang An (t - t n ) = Z &«0A'„(* - C). (19) 



ASYNCHRONOUS FM SPECTRA 2155 

The process on the left is the artificial process with independent a„'s and 
A„'s, whereas the process on the right hasthe6 n 's correlated, as observa- 
tion would require, and a different set of /„ 's representing the real process. 
Clearly the two processes are equivalent if one can find a transformation 
from the primed set of variables on the right to the unprimed set on the 
left. 

We observe the following characteristics of the two representations. 
The representation on the right demands that a transition occur at every 
/ = /„ for all n. To accomplish this b u must be different from b„-i , thus 
restricting the choice of the />„'s. The representation on the left admits 
independent choices of the a n 's, thus giving rise to masking of some 
transitions, since in fact, if a„-i = a„ there cannot be a transition at /„ . 
Furthermore, the set of /„ \s is a proper subset of the set of / n 's. To make 
the representation on the left useful, we must find how the two parameter 
sets transform. Toward this end define the following set of random 
variables 

X n =f(a n ,a n - l ) = l 1 ' a " " a '- i (20) 

(0, a„ ^ a„-i 

lorn = 1, 2, 3 ••• . 

Let Pj be the probability of obtaining a sequence of exactly jX's out. 
of N -\- j taking on the value unity. Then the probability P(N) of 
obtaining exactly N real transitions in T seconds in the process on the 
left of (19) is 

P(N) = Z/(A r +j)Py, (21) 

where f(N + j) is the probability of exactly N + ,;' transitions in the 
process on the left of (19). Equation (21) is a linear summand equation 
from which we would like to find a suitable /(• ) from the knowledge of 
P(-) and Pj. 

Not intending to make a general study of solutions of (21 ), one simple 
solution is presented in the next section for the case of exponentially 
distributed intervals A. In preparation for this discussion, we point out 
that if a„ is a discrete multilevel random variable with equally likely 
probabilities the set of random variables \X„\ in (20) is independent and 
therefore Pj is the binomial probability distribution. To demonstrate 
that indeed the set \X„\ is independent we have to show that the condi- 
tional probability distribution of X n given X n _i does not depend onX„_i . 
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With this in mind consider the joint characteristic function C(a>i , a> 2 ) of 
X„ and X n +i '• 

C(wi ,w 2 ) = (exp (io)iX n + iu2X n+ i))x„.x n + { 

= (exp [iuif(fln , a H -\)\ • exp \io) 2 f(a„ , a„+i )})«,„ .«„_,.«,,+, . 

From (20), fixing a„ and averaging first over a„+i and then over o„_i , we 
obtain 

C'(o), , ut) = ((exp [iuifian , a„-i )!)«.„_, 

•<exp [iwjf (an , a n+ i)}> a „ +l >a n ■ 

Since P{a n = y k ) - (1/Af ), k = 1, 2, • • • , M, independent of k we can 
write the last equation as 

k = M 

C(o), ,o) 2 ) = (1/Af) £ (exp {ictfi/(//*,a,,_i )}>„„_, , v 

fe-i (24) 

• (exp \iutf(Vk I fln+l)} >a ntI • 

Now 

(exp [iuifiyt , fl»-i) !)«,., = ^ exp (&*) + ( 1 - ^j , (25) 
and likewise 

(exp [iatfiyt , On+i) }>«.„+, = ^j exp (iw 2 ) + ( 1 - jj J ■ 

Since neither of the above averages under the summation sign in (24) 
depend on h; the joint characteristic function C(a>i , u 2 ) = C(wi)C(a> 2 ) 
which says that the random variables X n and X n -i are independent 
for all n. Clearly the above arguments still hold if a„ is allowed to 
take on a continuum of values. 

III. POISSON TRANSITIONS 

As a special case we assume that the number JV of t points in a fixed 
interval T obeys the Poisson probability law; consequently the prob- 
ability density of the intervals A between transitions is exponentially 
distributed, namely 

™- t Vi A<0. ™ 
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The characteristic function of A is then 

Ca(X) = (expr'XA) = [v/{v - ik)}. 
In particular it follows from (27) that 

vC A (\) - 1 



X = 



Ca(X) 



(27) 



(28) 



When (28) is substituted into (18) we obtain a very simple result for 
the spectrum, namely 



G+(w) = —Re 



y 



(C A (X)>> 



<C A (X)>x_ ' 



(29) 



where we made use of the fact that 

C A (X) 



i - Ca(x; 



x \X/ x 



which is purely imaginary. 

For black and white transmission a„ = ±1 with equal probability, 
and the spectrum reduces to 



A 2 

(?+(«) - — Re 



+ ^ 



2 j/ — z"Xi 2 v — i\ 



i-i 



2 i* - z'Xi 



1 
2 v 



iX 



2_ 



(30) 



where Xi = wj — w + a> c and X 2 = — m — u + co c from (12). By algebraic 
manipulation (30) is reduced to 

2 



0+(») " l 



V LOd 



03c)- 



(31) 



,/ [« d 2 - (« - W C ) 2 ] 2 + V 2 (C 

We see from this expression that the spectrum falls off as the fourth 
power of frequency. In a forthcoming section we shall present graphs of 
the various spectra. 

It is instructive to examine the physical meaning of the parameter v. 
This parameter is the average number of transitions per unit time of the 
virtual process. In fact the average number of transitions in the real 
process is i»[l — (1/M)], with the Poisson form of the density being 
preserved. To show that this is so, we observe that a solution of (21) is 
a Poisson probability distribution. In general, for M levels with Pr 
[ a „ = k] = l/M, k = 1, 2, • • • , M we have, from the previous section, 



I'j - 



IN_ 
J 



j\N\ \MJ \ M) 
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If f(N + j) is assumed to be 

-p N+j 

M+»=WTWV (32) 

we find from (21) that 

P(N) = E f(N + j)Pj = "^ (33) 

j=0 l\ '• 

where v\ = v[l — (1/M)]. Thus Poisson transitions with parameter v 
in the simple representation correspond to Poisson transitions with 
parameter v\ in the real process. 

IV. GRAPHICAL REPRESENTATION* 

In this section we present graphical results for the case of Poisson dis- 
tributed transitions. It is important to bear in mind the distinction be- 
tween the parameters for the virtual transitions, with which the calcula- 
tions are done, and the parameters of the real process, for which the 
results are reported. Here, we shall reverse the convention of (19) and use 
primes to distinguish "virtual" parameters. Frequency and fequency 
deviation are normalized to the average transition rate, i.e., 

/ _ M - 1 M -- 1 a - We 



M " M 2ttv 

(34) 

K' - M — l K - — - 1 — 
M M tv' 

The normalized characteristic function, with X„ defined as in (12), may 
then be written as 

1 



Ca(«0 = i o • taK' n ,\ ' (35) 

1 — Ziri 



^(-f-0" 



As our first example we consider a to be a discrete random variable 
taking on the possible values 2n — (M + 1), n = 1, 2, • • • , M with 
equal probability (1/M). Using these facts 

(C A (a)) a = x (1 - iZvp') + iirk'y, (36) 

where 



1 n=A 



M Z?i , , , n (aJC n ,V ' 



I + 4tt 



(<?-,)' 



ASYNCHRONOUS FM SPECTRA 

i r" 



2159 



V = m S 

.1/ »=1 



a» 



1 +4 



-pf-*f 



where we let a„ = In — (M + 1). Using (36) in (29) we obtain, for 
the virtual process 

^ + (/3') 1 - -v 



A 2 (1 - .f) 2 + 7r 2 (2.f/5' - K'yy 

and for the real process 

vG + (fi) M - 1 vG + {fi') 



^xi-1. 



[37) 



(38) 



A 2 .1/ A 2 

We have plotted this normalized spectral density (38) as a function 




1.0 1.5 2.0 2.5 

NORMALIZED FREQUENCY 



4.0 



Fig. 1 — Spectra for 2-level FM FAX .signals. 
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of the normalized frequency /3 for several values of normalized frequency 
deviation. Only the positive spectrum is shown since it is symmetrical 
about the normalized carrier, # = 0. 

The binary case is shown on Fig. 1. We note that these spectra con- 
tain none of the spectral lines which appear with synchronous modula- 
tion. There is, nevertheless, a tendency for the spectrum to be concen- 
trated about the frequency /3 = fc/2. Umiormalized, this is cj = w c ± 
Ktv. 

Higher level cases, M = 4 and 8, are shown in Figs. 2 and 3, respec- 
tively. These are similar to the binary case except that they have M 
levels and therefore M frequencies where concentration tends to occur. 
These frequencies are approximately /? = (2n — l)K/2, n = 1, 2, ■ • • , 
M/2. 
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Fig. 2 — Spectra for 4-level FM FAX signals. 
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Fig. 3 — Spectra for 8-level FM FAX signals. 

As a second example, we shall retain the Poisson distribution of transi- 
tions times, but allow the amplitudes to be continuously distributed over 
the interval [— r,r]. This is not true analog representation, but corre- 
sponds to "sample-and-hold" operation with exponential holding times. 
For this case the probability density of a is 

P(a) = l/2r, -r ^ a ^ r (39) 

and the expected value of the characteristic function becomes 



< c >°=£/' 



da 



v + id) — iiiidCb 



v , v + ibi — iu>dT 
= l n — In — r~- — r~- — 
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Fig. 4 — Spectra for continuous distribution FM FAX signals. 

Inserting this into (29) we obtain 

<? + (/3) R 2 + U + 2ttK)<po 



A"- R* + (** + 2vK)* ' 



(41) 



where 



R = In 



|~ 1 + (2tt/3 - 7rK) 2 T 
,1 + (2tt/3 + tK)*] 



<pa = arctan (2x0 — -wK) — arctan (2x/3 + tK). 

The normalized frequency deviation K is now modified to include r, 

namely 



K = rcod/irv. 



(42 
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For this continuous case, spectra for various K are shown on Fig. 4. 
The shape of these is very nearly rectangular, with height of \/2K 
and width of K/2, for the displayed positive spectrum. 

Considering the above results, together with results from a previous 
study on digital FM, 1 it is interesting to observe that in all the plots 
the shape of the spectrum is approximately the same as the first order 
probability density function of the baseband modulation process when 
K is large, in accordance with the adiabatic theorem. 2 
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